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Abstract. The possible control of competitive invasion by infection of the invader and multiplica¬ 
tive noise is studied. The basic model is the Lotka-Volterra competition system with emergent 
carrying capacities. Several stationary solutions of the non-infected and infected system are iden¬ 
tified as well as parameter ranges of bistability. The latter are used for the numerical study of 
invasion phenomena. The diffusivities, the infection but in particular the white and colored multi¬ 
plicative noise are the control parameters. It is shown that not only competition, possible infection 
and mobilities are important drivers of the invasive dynamics but also the noise and especially its 
color and the functional response of populations to the emergence of noise. 
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1. Introduction 

The main aim of modeling biological population dynamics is to improve the understanding of the 
functioning of food chains and webs as well as their dependence on internal and external condi¬ 
tions. Hence, mathematical models of biological population dynamics have not only to account for 
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growth and interactions but also for spatiotemporal processes like random or directed and joint or 
relative motion of species, as well as the heterogeneity of the environment. Early attempts began 
with statistics, exponential growth, physico-chemical (neutral) diffusion, and Lotka-Volterra type 
interactions. These approaches have been continuously refined to more realistic descriptions of the 
development of natural populations. 

Ecological and epidemiological models are known since more than 200 years. First attempts to 
merge these models appeared only about 30 years ago, cf. lUl Cl 121 [l3]] as well as ll55l l56ll . In¬ 
fectious diseases are prominent examples of biological invasions and continue to (re-)emerge in 
modern times. The negative econo-ecological effects of bioinvasions (HHI]] have led to a remark¬ 
able hype of bioinvasion research incl. modeling, cf. lfT5llJ7ll4^l^l59l . The history of research 
on stochastic processes and integration is long as well, historical surveys have been published, cf. 
ifTSllMll . The seminal work by Ito ifTTl and Stratonovich iBTIl should be particularly recognized. 
The present, to a large extent numerical study combines aspects of spatial eco-epidemiology and 
environmental stochasticity, namely the diffusive invasion of an alien species, its competition with 
the indigenous resident, and its biocontrol through targeted infection in a noisy environment. Con¬ 
trary to previous publications Il25ll26l . the environmental variability is modeled as external multi¬ 
plicative noise, in some cases with a certain functional response of the populations. 

Modeling environmental variability with multiplicative white noise goes back to the 1970s. 
Not only did May ll^ IbTII introduce the model that is used until today as a perturbation of the 
growth rate of a population by “white noise” but only a few years later, a more mechanistic basis 
of this model was developed. Branching processes provide a stochastic model that describes the 
number of offspring for a given number of individuals Zi within one generation i. The popula¬ 
tion number Zj+i in generation i -I-1 is updated for a given population Zi according to previously 
chosen probability distributions. This model was extended by Smith and Wilkinson ll421 l50ll by 
a stochastic process Q which modulates for each new generation the offspring that is generated. 
Branching processes in random environments (BPRE) provide an individual-based model for pop¬ 
ulation growth. For large population numbers, a BPRE can be approximated by a stochastic differ¬ 
ential equation (SDE) that accounts both for demographic as well as environmental stochasticity. 
Keiding ll20ll conjectured the form of the resulting diffusion approximation, his conjecture was 
rigorously proven by Kurtz ll2^ . The model introduced by May is obtained from Kurtz’ diffusion 
approximation by neglecting the term due to demographic stochasticity. Thus, the SDE model for 
environmental stochasticity is derived from the influence of a random environment on the popula¬ 
tion dynamics of a branching process but neglects its demographic stochasticity. 

In this study we consider the properties of environmental variability in more detail. The mul¬ 
tiplicative noise term implies that the effect of environmental fluctuations on the individuals of a 
population is additive. Whereas this seems reasonable for small population densities we propose 
that for large population numbers the effect of individual responses to environmental fluctuations 
on the population should decrease. Thus, we propose that the population-dependent response to 
environmental noise saturates for large population numbers similar to the functional response of 
predators at large prey densities and therefore we model the population-dependent response to 
environmental noise in a completely analogous way. 

The subject of our study is the influence of environmental stochasticity on a biological inva- 
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sion. In order to account for spatial spread we extend our system of SDEs by diffusion terms so 
that we obtain a system of stochastic reaction-diffusion equations. Spatiotemporal environmental 
fluctuations are represented by time-dependent random flelds. In contrast to previous studies, we 
consider random fields that are correlated both in space and in time, i.e., spatiotemporally colored 
noise. 

The assumption of uncorrelated white noise is usually justified by the coarseness of temporal 
or spatial scale, respectively. If spatial or temporal correlation length are much shorter than the 
time or length scale of interest, it seems valid to consider the time-dependent random field as 
uncorrelated. However, particular care must be taken when considering spatiotemporal dynamics 
driven by noise. Stochastic differential equations driven by uncorrelated noise can usually solved 
over a function space such as and this remains true if the system is extended to a reaction- 
diffusion equation over one-dimensional space (d = 1). But for spatial dimensions d > 2 solutions 
for stochastic reaction-diffusion equations driven by uncorrelated noise can only be guaranteed in 
a space of generalised functions, see e.g. (381 El]]. The reason for this phenomenon is that the 
Laplacian cannot smooth uncorrelated noise sufficiently for spatial dimensions exceeding 1 so that 
a solution may contain peaks resembling the S distribution. Not only is the physical significance 
of these solutions debatable but also numerical approximations are not capable of capturing this 
aspect of the continuous system. Here we take a pragmatic point of view on this difficult problem 
and present numerical solutions for temporally and spatially white noise where the space may be 
interpreted as a discrete lattice whose nodes interact by a discrete Laplacian. 

2. Resident-invader competition with infection in the invader 
population 

For the invasion of a resident population by a competing invader, the Lotka-Volterra competition 
model is used, i.e.. 


= riA^i - ciiA^f - C12N1N2, 


( 2 . 1 ) 



( 2 . 2 ) 


where Ni and N 2 are resident and invader respectively. Carrying capacities will not explicitely 
be introduced because they can suppress a higher variety of solutions and rather appear as special 
cases ]l8l|22l|28l|46|]. The r’s stand for the growth rates and the c’s for the inter- and intraspecific 
competition. 

A specific infection of the invading population can be used as biocontrol measure to stop and 
reverse the invasion, cf. (SJImKT^ISS]]. To model this, the invader population is split into suscep- 
tibles S and infecteds /, 


A^2 = ^ + /. 

Then, the local dynamics reads with notation X = {Xi = Ni, X 2 = S, X 3 = 1} 
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(2.3) 


-^ = fi (X) = riXi - c^,X! - c, 2 X^{X 2 + X 3 ), 

dX X X 

-T^ = f 2 (X) = raXa - C 22 X 2 (X 2 + X 3 ) - C 21 X 1 X 2 - A (2.4) 

at [Ji 2 + ^3 

dX X X 

-X = h (X) = r2X3 - C22X3(X2 + X 3 ) - C 21 X 1 X 3 + A .^ ^ - ijXs , (2.5) 

at [A 2 + A 3 

where A is the transmission coefficient of the disease and jj, the disease-induced higher mortality 
rate of the infecteds. The exponent q allows to describe mass-action type (q = 0) and frequency- 
dependent transmission (q = 1) of the disease respectively ifTlI^. 


However, one cannot expect that growth rates and competition intensities of susceptibles and 
infecteds are the same. They should rather be split and could be ordered like 

r2 {r2,r3} & ^3 < r2, 

C 12 ^ {ci 2 ,Ci 3 } & Ci 3 < C 12 , (2 6) 

C2I ^ {c21,C3i} & C21 < C31 , 

C22 ^ {C22, C23, C32, C33} & C23 < C33 < C22 < C32 . 

The ordering of the intra- and interspecific competition coefficients of susceptibles and infect¬ 
eds depends on the biological species, cf. However, it can be certainly accepted ll46ll that 


C 23 ^ C 22 A C 33 < C 32 . 


System ( |2.3| - |2.5| ) then changes to 

dX 

= /i(X) = (n - cnXi)Xi - (C 12 X 2 + ci3X3)Xi, (2.7) 

dX X X 

= /2(X) = (r2 — 022 X 2)^2 — (C 21 X 1 -\- C2zXz)X2 — X ^ ^ '^q-> 

dX X X 

= fs(X) = (n-l^- C33X3)X3 - (C31X1 + C 32 X 2 )X 3 + A ^ ^ . ( 2 . 9 ) 

at (A 2 + A 3 p 


For convenience, the model of the local dynamics is not analysed in terms of Xi, X 2 and X 3 
but rather in Xi, i and N 2 = X 2 + X 3 where i is the prevalence, i.e., the infected fraction of the 
total invader population Xaffll, 

. ^ X 2 ^ X 2 
X 2 + X 3 “ N 2 


Having in mind that 

di 1 / dX2 .dN2\ 

dt N 2 \ dt ^ dt ) ' 


( 2 . 10 ) 
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it follows 


^ = riXi - cnX2 - [ci2(l - i) + cigi] , (2.11) 

di 

^ = {^3 ~ + -^-^2 ^ ~ /^ + -^2 [(C 22 ~ C32)(l — i) + (C 23 — £ 33 )^] + Xi [C 2 l — C 3 i]| * 

*i{l-i), ( 2 . 12 ) 

^ = {G'(7V2, i) - [C2i(l - i) + C3iz] Xi} X 2 , (2.13) 

with 


G{N2,i) = r2(l - i) + (r3 - fi)i - N 2 [ 032(1 - if + (023 + C32)i(l - i) + . (2.14) 

The latter expression is also found for predator-prey systems with infected prey Il4^ . Note that 
if the resident Xi resp. the predator in ll46ll cannot distinguish between susceptible and infected 
invader resp. prey, the temporal change of the prevalence becomes independent of the type of 
interspecific ecological interactions such as competition and predation. It only contains terms 
describing the intraspecific competition of susceptibles and infecteds in the infected population. 

2.1. Stationary solutions and stability for frequency-dependent (standard) 
incidence q=l 

In phytopathology, the transmission of especially fungal diseases is described with standard inci¬ 
dence Il5^ . A corresponding model of the invasion of a fungal disease over a vineyard has been 
investigated in BH]. Further on, only the standard incidence is considered, i.e., q = 1. 

The infection-free system, i.e., i = 0, N 2 — X 2 , is the Lotka-Volterra competition model with 
its known stationary solutions and their stability ranges. Especially interesting for the consideration 
of spatial invasions is the bistable parameter range 


riC22 - r2Ci2 > 0 A raCn - riC2i > 0, 


Ti T2 

when both the invader-free (—, 0 , 0 ) and the resident-free ( 0 , 0 , —) states are stable and can 

Cll C22 

compete for space. The opposite case may also happen: The invader arrives already infected, i.e., 
z = 1 , X 2 = X 3 , and the invader-free (—, 0 , 0 ) and the resident-free ( 0 , 0 ,-) states can be 

Cll C 33 

both at once stable for 


?’lC33 - (?’3 - /i) Ci 3 > 0 A (r 3 - fl) Cll - riC 3 i > 0 . 

Ti 

However, the latter as well as the possible bistability of resident-only (—, 0,0) and (susceptible- 

Cll 

infected)-invader-only ( 0 , states will not be considered here. 
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2.2. Spatiotemporal dynamics in a variable environment 


The main focus of this study is to consider the spatiotemporal effects of a more detailed model 
of environmental variability. We assume that all species spread randomly so that mobility can be 
described as diffusion with coefficients D = {Du = Di ; Dij = 0 Vi 7^ j ; i, j = 1,2,3}. Also we 
add Gaussian random fields ^(r, t) = f); i = 1,2,3} to system (2.7-2.9) so that we obtain 
the system of stochastic partial differential equations 


- DAX(r, t) = f [X(r, t)] + g [X(f, t)] ${f, t ), (2.15) 

where the matrix function g(X) = {gij{X.); i,j = 1, 2,3} determines the density-dependent 
noise intensity. We consider horizontal processes with position vector f = {x,y} and corre¬ 
sponding Laplace operator A = 8 “^/dx^ + 8 “^/dy^. In literature, often temporally and spatially 
uncorrelated “white” Gaussian fields with zero mean and delta correlation have been considered 


=0, ( 6 (ri,fi)^*(f 2 ,t 2 )) = S{fi - f2)S{ti - ^ 2 ), i = 1,2,3. (2.16) 

Here, we investigate the effect of extending this model by correlated “colored” noise with 
correlation lenghts r and A in the temporal and spatial domain, respectively. Apart from using 
colored noise we also investigate a generalisation of the density-dependent noise g (X). Purely 
diagonal, linear multiplicative noise 

guiX.) = uJiiXi ; gij{X) = 0 V z 7 ^ j ; z, j = 1, 2, 3 ; (2.17) 

can be interpreted as a model where individuals respond independently to stochastic environ¬ 
mental variability. Thus, the effect of environmental fluctuations on each individual directly trans¬ 
lates into variability at the population level - the response at the population level is additive. The 
alternative model suggested here is based on the assumption that in large populations individu¬ 
als do not repond to fluctuations independently from each other. Instead we propose that larger 
populations respond to environmental variability in a more robust way, i.e., neither favourable nor 
adverse effects influence the population proportional to the number of individuals: 

(jj. . 

gi,{X) = --; z ,7 = 1,2,3; 1 < m < rz < 2 . (2.18) 

^ O'ikXjl 
k=l 

For m = n, the parameter Uij/aij is the maximum noise intensity that is reached asymp¬ 
totically for large populations Xj. The parameter 7 *^ is the population level at which the noise 
intensity reaches half of the maximum level coijjaij. Thus, this parameter describes the ability of 
a population to collectively reduce the effect of noise - the higher 7 ^^, the higher the population 
level must be until population Xj is appreciably affected by environmental variability. 

For m < n, the noise intensity even decreases and eventually vanishes for high population 
densities. However, these values are never reached. 
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In previous papers (25112^, it was shown that a certain variability of the environment and the 
mobilities of the competitors are the system-driving forces. Extreme events such as landslides lead 
to bare ground re-invadable by both resident and alien species. These events at random times, size 
and locations are not considered here. They are replaced by white and colored noises |l42l|47l|48]]. 
Again, the biocontrol of the invasion through a specific infection of the non-indigenous species is 
studied. 


3. Numerical methods 


The numerical solution of stochastic partial differential equations is a difficult problem and the 
subject of current research. For this reason we explain how the spatiotemporal model ( |2.15| ) can 
be solved numerically and how spatiotemporally correlated noise can be generated. We follow a 
finite difference approach where in a first step the spatial domain is discretised. In this way the sys¬ 
tem of stochastic partial differential equation is approximated by uncoupled stochastic differential 
equations that are solved numerically in a second step. For the first step, we use the semi-implicit 
Peaceman-Rachford method |l40l |52l which, in particular for stochastic equations, often seems to 
be more robust than the simplest explicit scheme. In the following we will explain how the dis¬ 
cretised system can be solved using the derivative-free Milstein method and how spatiotemporally 
correlated noise can be generated. 


3.1. Derivative-free Milstein method 

For numerical integration, the derivative-free Milstein method is used, cf. 112711^1^ but also the 
short descriptions |[l0l|44]|. Sometimes and in particular for the purpose of this study, it is even 
sufficient to consider purely diagonal intensity matrices 


Pn(X) 




; ajj = l;i = l,2,3; l<m<n<2. 


Then, the Milstein scheme reads with time step At and Stratonovich interpretation 


^;+A. = x; + /.(A1)A( + ft.(A7)AWj + ^ [ft.(A7) - fti(A;] (AWj)", 

with 

Xi = Xi + MXl)At + guixi)^t, 

and 

AVF, = WUt - Wi - y/AtXiO, 1). 


(3.1) 


(3.2) 


As usual, AA(0,1) stands for the normal distribution with zero mean and unity variance. The 
required uniformly distributed random numbers are generated with the Mersenne Twister ll29l . the 
normally distributed with the common Box-Muller algorithm Q. 
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3.2. Generation of correlated Gaussian random fields 


Garcfa-Ojalvo and Sancho ifTTl [T2l developed a method for generating spatially and temporally 
colored noise from the stochastic reaction-diffusion equation 


dC,{f,t) . 1 X 

—— = — AC(r,t) - + -r){r,t) 

at r r T 


(3.3) 


Here, the term r^(r, t) stands for uncorrelated (white) noise. The parameters r and A determine 
the correlation lengths in the temporal and the spatial domain, respectively. In addition, a scaling 
factor e for the variance of normally-distributed random variables in Fourier space has to be cho¬ 
sen. The authors introduce Eq. (|3.3|) as an analogon of the Ohmstein-Uhlenbeck process which 

is the solution of ( |3.3| ) without the spatial term —But it has to be noted that for two- 

dimensional space (and spatial dimensions exceeding two) the solutions of equation p.3| ) will be 
generalised functions rather than functions in a space such as Thus, strictly speaking, we are 


generating discrete spatiotemporally colored random fields that are derived from the model ( |3.3| ) 
without being approximate solutions of the continuous problem. 

The spatiotemporal random field is simulated by transforming a discretised version Qj (t) 

of ( |3.3| ) to Fourier space: 


CuAt) = {Axf exp(-ifc • r)Oi(^) 


(3.4) 


We denote the discrete Fourier transform with greek indices /U, v rather than i, j and the 

coordinate k in frequency space is 




//, Z2 = 0,..., L — 1. 


Then the Fourier transformed field at the next time step t + At is calculated by 


+ At) = Cf,^{t) exp (— ^At ) -b 


' e{LAx) 


2 r 


TG 




1 — exp ( — 


a 


jiu 


(3.5) 


Here, is the Fourier transform of an uncorrelated Gaussian noise field aij. The efficiency 
of this method is increased by directly generating the Fourier transformed field The (complex¬ 
valued) discrete Fourier transform of real-valued fields a^j obeys some symmetries that lead 
to the following restrictions: 


0/-IXU — ^L-n,L-v 


OLfxi! G 


for ii,u = 0,- 


(3.6) 

(3.7) 


where denotes the complex conjugate of e C. The condition (3.6) means that have 
the same real part as found by reflecting through the centre {[i = p = 1/2) whereas the 

imaginary parts only differ by opposite signs. 
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Also in p.5| ), is the Fourier transform of the discretisation of the differential operator 
id-A^A: 

2A^ r / 27r/u\ ( 27rz/' 


^[JLV 1 


(Ax)^ 


cos 


+ cos 


For further details on the derivation of this method, we refer the reader to ifTTll or ifT^ . 


L = 
(3.8) 


4. Resident-invader competition-diffusion model with infected 
invader in a variable environment 

4.1. Local dynamics with multiplicative white noise and induced transitions 

Not surprising and like for g = 0 (531153, Lehmann (^ and Woyzichovski (6011 found disease- 
induced oscillations fox q = 1 as well. Their interesting result was that there may be bistability 
of the resident-only state and the limit cycle when coexisting resident, susceptible and infected 
invaders oscillate, cf. Figure[^ The following parameters have been used for the latter setting: 

ri = 1.2500, ra = 1.0000, rg = 0.6775 

cn = 0.5000 , ci 2 = 0.8500 , Cai = 0.4250 , Cis = 0.5000 , C 31 = 0.5000 (4.1) 

C 22 = 0.6000 , C 23 = 0.4000 , C 32 = 0.6010 , C 33 = 0.5000 
A = 0.4000, M = 0.0250 

All other semi-trivial states turn out to be unstable for this parameter range. 




Figure 1: Bistability of resident-only state and oscillating coexistence 


Now, one effect of white noise in locally multiple stable systems is shown, i.e., the switch from 
one stable attractor to the other for sufficiently but not too high noise intensity. For simplicity, the 
linear density dependence ( |2.17| ) of the intensity is chosen, which has been successfully applied to 
numerous cases. Here, the leaving of the initial limit cycle is demonstrated. 


The three subfigures of Figurej^ show typical outcomes of hundreds of simulations with dif¬ 
ferent seeds of the random number generator. The left subfigure shows the persistence of the 
limit cycle whereas the middle demonstrates the expected leaving of the cycle for the other stable 
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Figure 2: Noise-induced transitions from oscillating coexistence to resident-only state resp. resi¬ 
dent extinction. Linear density dependence ( |2.17| ) with a;jj = 0.1 i = l,2,3. 

stationary solution, i.e., the resident-only state. The result in the right subfigure appeared a bit 
unexpected, however, due to some catastrophic shift, the resident died out and the remaining sus¬ 
ceptible and infected invader survived. One should have in mind that the latter (0, X|) state is 

unstable to the reintroduction of the resident. 


4.2. Invasions and noise I 

4.2.1. Linear noise and biocontrol of invasion 

For the beginning, the results in the mentioned previous papers ll^l26l . where simulated landslides 
led to bare land competitively re-invadable by resident and invader, are reproduced with external 
noise. The initial condition is a “red” invader patch at its emergent carrying capacity r 22 /c 22 at the 
“upper left comer” of the “green” habitat of the native species at rn/cn. This patch should exceed 
the related critical patch size. Otherwise it will simply decay regardless of its competitive strength 
and mobility llTTlI^ . Zero-flux boundary conditions are applied. 


The linear noise (2.17) and parameters 


from the previous publications have been taken: 


ri = 1.000, r2 = 1.000, — 0.800 

Cii = 1.000, Ci2 = 1.300, C 21 = 1.200, C 13 = 1.299, 
C 22 = 0.999 , C23 = 0.998 , C 32 = 1.001, 
A = 0.405 , = 0.200 

Uii = 0.250 ; A = 45.000 ; i = 1,2,3. 


C31 = 1.201 

C 33 = 1.000 


(4.2) 


The invader patch spreads and seems to grow unstoppable. 

Then, a biological control measure is applied. The invader population is partly infected and the 
invasion successfully rolled back. 

The changes of the fraction of the invaded area can be seen in Figurej^ 
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Figure 3: Simulation 1: Parameters as given in ( |4.2| ). The first row shows results with white noise, 
the second row shows spatiotemporally colored noise with uju = 0.03, e = 0.001 and correlation 
lengths T = 1 and A = 1. 



t=1200 1400 1600 1800 2000 

Figure 4: Simulation 1 continued: Partial “blue” infection of the invader at t=1200. 


4.2.2. Nonlinear response to noise and noise control of invasion 


We now consider the saturating response to noise described above ( |2.18| ) with 

TTL = n = ‘2 ] iOii = 50 , ^22 = 0.1 , CJ 33 = 10 5 
7ii = 100, 722 = 3 , 733 = 300. 


(4.3) 


The extinction of the invader due to hostile environmental conditions is shown in Figure]^ The 
resident is used and adapted to the environment and happily survives. 


This dynamics is only due to the specific noise response of the populations. All growth and 


interaction parameters remained the same as in sec. 4.2.1. 


11 









Figure 5: Fraction of invaded area before and after partial infection. 



t=100 300 400 500 700 



t=500 1000 1500 2000 2500 


Figure 6 : Repressed invasion through population-specific nonlinearly density-dependent noise in¬ 
tensities for white noise in the first row. In the second row for colored noise with correlation 
lengths r = 1 and A = 1, cjn = 0.05, U 22 = <^33 = 0.01 qualitatively similar behaviour is ob¬ 
served, albeit on a much slower timescale. However, other simulations have shown that stronger 
temporal and spatial correlations can support invasions. 

4.3. Invasions and noise Ila 

Coming back to the parameter range of bistability of resident-only state and oscillating coexistence 
of all three populations, i.e., parameters ( |4.1| ), the initial condition is chosen as uniformly popu¬ 
lated by the resident Xi at its emergent carrying capacity. At a defined location at the boundary, an 
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initial patch of the invading populations X 2 and attempts to spread. Again, simply linear noise 
( |2.17[ ) is applied and the influence of increasing noise intensity ccjj = cc ; i = 1,2,3; studied. 

It is seen in Figure]^ that a successful invasion requires a certain supercritical noise intensity. 


(D 

■o 

CtS 

> 



Figure 7: Decline resp. growth of invaded area for increasing white noise intensity. 

Somehow, the resident supports the invasion of its own area. Diffusion and noise enhance the 
mixing of resident and invaders at the front. Therefore, all three together jump on the stable limit 
cycle of coexistence and invade the remaining invader-free area, cf. Figure 
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Figure 8 : Upper row: Successful fast invasion of the resident’s area for linear white noise ( |2.17| ) 
with Wii = 0.5 ; i = 1, 2, 3. Lower row: Stopped invasion of the resident’s area for nonlinear noise 
(|43]) with u)ii = 0.25 ; i = 1 , 2 , 3 ; 711 = 1.00 , 722 = 25.00 , 733 = 4.00. 


The following Figurej^ shows that weakly correlated noise still allows for invasion but stronger 
correlated does not. 



200 500 2000 4000 



t=0 too 150 500 1200 

Figure 9: Same parameters as in Fig.[^ The invasion is successful for colored noise and correlation 
lengths T = A = 1 but unsuccessful for r = A = 15. Other parameters e = 0.001, cou = 0.05. 

After some difficulties at the beginning that can be the end for the invader at lower noise 
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intensities, the purple invader patch turns into the blue of the limit cycle. Finally, the resident 
survives but has to share its habitat with the aliens. However, if the nonlinear response or colored 
noise is applied, the invasion can be stopped and rolled back again. 

4.4. Invasions and noise Ilb 

Now, it is assumed that the limit cycle of resident and invaders has already invaded most of the 
area and only a small part is left for the resident alone. Again, the parameters ( |4.1| ) are used. For 
this situation several interesting patterns appear that are again purely due to different properties of 
the environmental noise. 

4.4.1. Dependence on noise intensity 

It is surprising that the resident turns out to be strong enough to defeat the invasion as far the noise 
is below a subcritical threshold. The fraction of invaded area over time is plotted in Figure[T^ An 



Figure 10: Decline resp. growth of invaded area for increasing white noise intensity, 
example for the defeat of invasion is given in FigurefTTj 
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t=0 1500 3000 4500 6000 


Figure 11: Finally successful defeat of invasion for no noise (upper row) resp. white noise with 

LUii = 0.1; i = 1,2,3. 

A preliminary conclusion is that increasing linearly density-dependent white noise supports 
invasion. In Figure [T^ the cloudy result for cou == 0.25 is shown. 



t=100 1000 2000 4500 6000 


Figure 12: Cloudy invasion for linear white noise with uu = 0.25 ; i = 1,2,3. Initial condition as 
in Figure[TT| 

Increasing the noise intensity and both correlation lengths (r = A = 20) leads to a situation 
where native and resident gain and lose control over parts of the spatial domain in an alternating 
fashion (Figure p^. It seems that the resident is slowly getting the upper hand: att = 7500 roughly 
two thirds of the domain are occupied by natives, at t = 10000 the resident has lost a few areas 
and displaced the invader in a few others. 
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t=50 500 2000 7500 10000 


Figure 13: Resident and invader fight under colored noise for territory creating a highly irregular 
pattern for a;** = 0 . 1 , e = 0 . 001 , r = 20 , A = 20 . Initial condition as in Figure[TT| 

4.4.2. Metapopulation patches 

This does not necessarily change for nonlinear noise, however, one setting is found where the 
fraction of invaded area is dropped down from initially 98% to 21%. The resident population splits 
into three metapopulations, spatially separated by the invader populations, cf. Figure[T4| 



t=500 1000 1500 2300 6000 


Figure 14: Stable formation of spatially separated patches of the resident population for white 
noise intensity uju = 0.25 ; z = 1,2,3; 711 = 1.0, 722 = 25.0 , 733 = 4.0. Initial condition as in 
FigurepTj 

As observed above, higher noise intensities oj help the invader to establish by stabilising the 
coexistence limit cycle. In contrast, stronger correlations i.e. increasing correlation lengths r in 
time or A in space generally enable the native species to displace the invaders. Only for relatively 
high correlation lengths of 10 or above the native species is able to form patches and avoid a cloudy 
invasion as in Figure [T^ 

In Figure [T^ for a low noise intensity tUjj = 0.03, r = 1 and A = 15 we observe the emergence 
of a quasi-stationary pattern similar to Figure [T4| 
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t=340 1000 2000 3000 6000 


Figure 15: Resident slowly displaces invader while forming seemingly stable spatially separated 
patches, similar to the results in Figurebut without nonlinearity in the noise term. Parameters: 
iOii = 0.03 , e = 0 . 001 , T = 1 , A = 15. Initial condition as in FigurefTTj 

4.4.3. Spiral waves 

Wavy structures are found as well, however, at the cost of full invasion. One example is presented 
in FigurefT^ 



t=500 1700 2250 5500 6000 


Figure 16: Stable formation of spatially separated patches of the resident population for linear 
white noise intensity = 0.3; i = 1,2,3; 711 = 0.5625, 722 = 4.0, 733 = 9.0. Initial 
condition as in Figure[TTj 

For nonlinear noise, the spiral waves seen in Figure [T^ can also be observed if the noise is 
colored. For small values of spatial and temporal correlation lengths (r = A = 1), smaller and 
slightly more irregular spiral waves can be observed, cf. Figure [T7| 

For increased correlation in time or space (e.g. r = 1 as above but A = 5), the native population 
can still maintain a spreading front that eventually displaces the spiral waves formed in its wake, 
cf. Figure [T^ 
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t=1000 t=2000 t=5000 t=7000 t=9000 


Figure 17: Formation of spiral waves - similar to Figure [T^ but less regular - for colored noise 
intensity uou = 0.03; i = 1,2,3; e = 0.001, r = 1, A = 1, 711 = 0.5625, 722 = 4.0, 733 = 
9.0. Initial condition as in Figure[TTj 



t=1000 t=2000 t=5000 t=7000 t=9000 


Figure 18: For increased spatial correlation the spiral waves are displaced by a front of natives. 
Parameters: ua = 0.03 ; z = 1, 2,3 ; e = 0.001, r = 1, A = 5, 711 = 0.5625 , 722 = 4.0 , 733 = 
9.0. Initial condition as in Figure[TTj 

5. Conclusions 

Populations are exposed to fluctuations of many environmental parameters such as nutrient avail¬ 
ability, temperature etc. that may have positive or adverse effects. Whereas it would be impractical 
to explicitly consider a host of factors that each on their own may only have a small influence on 
growth or decline of the population it is possible to represent the collective effect of these factors as 
stochastic environmental variability. The standard model for stochastic environmental variability 
in population dynamics are stochastic differential equations with a multiplicative noise term. 

In this model, the matrix of maximum noise intensities u is the only parameter that can be used 
for capturing all aspects of environmental fluctuations. The standard assumption that environmen¬ 
tal noise is temporally and spatially uncorrelated neglects the fact that many environmental factors 
are, in fact, typically correlated. 

Whereas at first glance it seems that this can be convincingly justified by assuming that the 
spatial and temporal correlation lengths r and A of the noise are much shorter than the spatiotem- 
poral scale under consideration, we have demonstrated in this study that using correlated instead 
of uncorrelated noise may lead to qualitatively very different model behaviour. Thus, neglecting 
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possible correlations may, in fact, lead to different explanations of the observed system behaviour. 

Another implicit assumption that has previously been unquestioned is the linear increase of 
noise intensity with population number. This model implies that environmental effects on each 
individual in a population simply add up to an overall effect on the population. In our opinion 
it is likely that for increasing population numbers, perturbations should not independently affect 
each individual but rather saturate due to interactions of the individuals so that the collective re¬ 
sponse of the population saturates to a maximum noise intensity for large population numbers. 
This model requires an additional matrix 7 which characterises the ability of the population to 
“buffer” stochastic fluctuations: for low values of '^u, population Xi is exposed to intensities close 
to the maximum noise level oju for low or moderate population numbers whereas for a population 
with a large 7 ** the noise intensity uju is only reached for high population numbers. 

The results presented in this paper have been obtained for purely diagonal noise intensity ma¬ 
trices. More complex forms are left to future work. 
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